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Abstract 

A  new  approach  to  guidance  algorithm  design  using  the  concept  of  the 
information  sets  (IS)  is  considered.  As  an  example  of  such  approach,  the  model 
of  defense  scenario  against  reentering  ballistic  missile  is  investigated.  In  this 
scenario,  a  maneuverable  decelerating  target  is  to  be  destroyed  by  a  hit-to-kill 
interceptor  using  an  IR  array  seeker  and  lateral  impulse  thrusters.  The  key 
element  of  the  proposed  approach  is  a  description  of  the  ’’interceptor-target” 
system  state  by  means  of  IS.  The  numerical  methods  for  fast  transformations 
of  the  IS  are  developed. 


2 


Contents 


1  Introduction  4 

2  Problem  statement  5 

3  Information  sets  construction  7 

4  Auxiliary  impulse  control  problem  9 

5  Solution  of  auxiliary  control  problem  11 

6  The  guidance  algorithm  against  a  maneuvering  target  based  on  the 

IS  concept  14 

6.1  Simple  one-dimensional  example .  14 

6.2  Problem  of  guidance  to  a  maneuvering  target  .  15 

7  Monte  Carlo  test  of  the  algorithm  of  guidance  to  a  non-maneuvering 

target  17 

7.1  A  priori  data  .  17 

7.2  Evaluation  of  the  guidance  algorithm  performance .  19 

7.3  Computational  properties  of  the  algorithm .  19 

7.4  Resume  .  27 

8  Monte  Carlo  test  of  the  algorithm  of  guidance  to  a  maneuvering 

target  28 

8.1  A  priori  data  .  28 

8.2  Evaluation  of  the  guidance  algorithm  performance . 30 

8.3  Computational  properties  of  the  algorithm .  31 

8.4  Results  of  the  Monte  Carlo  simulation  .  31 

8.5  Resume  .  32 

9  Conclusion  33 

10  Outlook  33 


3 


1  Introduction 


The  promising  approach  to  the  analysis  of  many  problems  of  control  with  incom¬ 
plete  information  has  been  developed  in  recent  years.  This  approach  is  concerned 
with  describing  positions  of  the  controlled  system  by  so-called  information  set  (IS), 
all  points  of  which  are  compatible  with  the  system’s  dynamics,  control  history,  mea¬ 
surement  model  and  observation  history  [1],  [2].  In  the  framework  of  this  approach 
the  original  control  problem  is  reduced  to  a  problem  of  control  of  the  IS  evolution. 
Usually,  the  aim  of  such  control  is  to  minimize  some  performance  index  as  a  function 
of  information  set. 

This  approach  is  especially  suitable  for  problems  of  guidance  with  nonlinear 
(bearings/range,  bearings-only)  measurements  corrupted  by  additive  non-Gaussian 
noise  distributed  within  given  boundaries.  It  is  well  known  that  in  such  situations 
a  Kalman-type  estimator  incorporated  in  the  guidance  loop  usually  fails  for  this 
measurement  model,  and  overall  guidance  performance  deteriorates  significantly. 

Using  proposed  IS-based  approach,  the  three-dimensional  problem  of  guidance 
with  incomplete  information  is  investigated.  The  impulse-controlled  interceptor  M 
is  guided  to  the  maneuverable  decelerating  target  T.  During  the  guidance  the  in¬ 
terceptor  M  performs  passive  angle-type  measurements  corrupted  by  additive  uni¬ 
formly  distributed  noise.  The  aim  of  the  guidance  is  a  minimization  of  the  relative 
distance  to  ensure  target  impact.  An  important  feature  of  the  problem  is  the  fact 
that  the  interceptor’s  control  impulses  are  orthogonal  to  its  longitudinal  axis,  so  the 
longitudinal  component  of  the  interceptor’s  motion  is  uncontrolled. 

Similar  problems  were  examined  in  works  [3],  [4]  in  the  context  of  minimax 
filtration  and  the  theory  of  differential  games. 

In  the  frame  of  this  contract  two  different  models  of  the  target’s  relational  motion 
are  considered: 

1.  lateral  acceleration  in  chosen  coordinate  system  is  constant  but  unknown; 

2.  lateral  acceleration  can  change  instantly  remaining  within  given  boundaries. 

In  both  cases  the  longitudinal  motion  of  interceptor  and  target  is  a  deceleration 
caused  by  the  drag  forces.  The  drag  coefficients  are  assumed  to  be  known. 

First,  we  consider  the  model  No.l  with  constant  lateral  acceleration  caused  by 
a  drag  force.  More  precisely  the  lateral  acceleration  is  approximately  constant.  It 
changes  much  slower  in  the  first  model  in  contrast  with  jump-wise  aerodynamical 
maneuvers  in  the  second  one.  After  formulation  of  a  mathematical  model  of  the 
guidance  problem,  we  define  the  information  sets  for  the  considered  dynamics.  To 
facilitate  numerical  operations  with  the  introduced  IS  we  propose  to  use  so-called 
Grid-Polygon  Approximation  of  the  IS. 

As  a  criterion  for  the  IS  evolution  control  the  miss  averaged  over  IS  is  taken.  The 
interceptor’s  control  is  chosen  in  the  class  of  so-called  ’’conditionally  programmed 
controls”,  i.e.  programmed  controls  which  are  updated  after  each  measurement. 
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Recently,  this  approach  has  been  applied  to  similar  guidance  problems  in  [5,  6] 
and  very  promising  results  have  been  obtained. 

To  update  the  control  the  auxiliary  optimal  impulse  control  problem  is  solved.  In 
this  problem  a  performance  index  has  the  meaning  of  expense  of  impulse  control  over 
time-to-go,  and  a  terminal  condition  means  zero  average  miss.  The  auxiliary  prob¬ 
lem  is  transformated  to  the  quadratic  mathematical  programming  problem  which  is 
solved  using  Kuhn- Tucker  optimality  conditions. 


2  Problem  statement 

The  guidance  geometry  is  shown  in  Fig.l. 


Figure  1:  The  guidance  geometry.  1  -  interceptor,  2  -  IR  array  seeker,  3  -  lateral 
impulse  thrusters 


It  is  assumed  that  the  target  T  moves  with  zero  attack  and  sliding  angles  and 
the  lift  force  can  be  neglected.  Then  the  target’s  motion  is  a  deceleration  caused  by 
the  drag  force: 

r-t  n 

l/i/i  =  Vo  -  or  /  V2(s)ds,  or  =  1  .  (1) 

Jto  IniT 

Here  Vr{t)  -  target’s  longitudinal  velocity,  -  target’s  mass,  p  -  atmospheric 

density,  S  -  target’s  frontal  area,  Cx  -  drag  coefficient.  The  coefficient  or  is  assumed 
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to  be  known.  The  relational  motion  of  the  interceptor  M  and  the  target  T  in  the 
reference  system  OXYZ  fixed  to  the  interceptor  is  described  by  the  equations: 


x(t) 

y(t) 
z(t ) 


/  to 


kx  Vt(s)  ds, 


to 


ky  Vt(s)  ds , 
kz  VT(s )  ds, 


v,[t)  =  E 

t ^  t  t  t 


Vm{1- )  =  Wq  —  aM 


to 


VM(s)  ds, 


(2) 

(3) 

(4) 

(5) 

(6) 


where  t0  -  the  initial  instant  of  time,  Vm,  Vy ,  Vz  -  interceptor’s  velocities,  (iq,  u2)  = 
(ui(t),u2(t))  -  interceptor’s  impulse  controls  represented  by  a  totality  of  impulses 
uik  =  ui(tk),  u2k  =  u2(tk)  which  are  applied  at  fixed  instants  tk  =  t0  +  kA,  k  = 
0,  1, .  .  and  A  -  given  time  interval.  The  constants  kx,  ky  and  kz  equal: 

kx  =  cos  tf  cos  tf),  ky  =  sin  kz  =  cos  <p  sin  ip. 


The  initial  parameters  x0,  y0,  z0,  Vo  and  coefficients  ky,  kz  are  assumed  to  be 
unknown  but  restricted  by  the  given  inequalities. 

During  the  guidance  the  interceptor  performs  measurements  (Cifc,  (2 k)  described 
by  the  equations 

C  ytpk)  ,  _  z(tk )  ,  _  7  n  1 

C lk  ~  ~7TT  +  hi k,  (,2k  —  ,  ,  +  1]2k ,  k  —  0,1,  ...  , 

x(tk)  X(tk) 

where  rj \k,  r/2k  -  uniformly  distributed  measurement  noises: 


Ihlfcl  <  C,  Ih2fc|  <  C, 


c  -  given  constant.  The  observation  process  and  control  process  are  terminated  at 
an  instant  tj,  where 

tj  =  inf{t  :  x(t)  <  e}.  (7) 

Here  e  -  given  small  distance,  e  >  0.  The  physical  meaning  of  the  condition  (7) 
is  that  target  observations  are  impossible  when  the  distance  ’’interceptor-target”  is 
lower  than  e  due  to  the  large  angular  size  of  the  target  and  overload  of  the  seeker’s 
IR  receiver. 


6 


3  Information  sets  construction 


By  definition,  the  information  set  (IS)  consists  of  all  points  of  the  phase  space 
whose  coordinates  are  compatible  (in  geometrical  sense)  with  dynamics  equations, 
measurement  model  and  restrictions  imposed  on  an  initial  state  vector. 

For  the  problem  under  consideration,  the  IS  I(t)  is  6-dimensional  one: 

I(t)  C  {x,VT,y,z,ky,kz}, 

assuming  the  coefficient  kx  to  be  known. 

Note  that  the  dynamics  equation  described  the  motions  along  Y  and  Z  axis  as 
well  as  the  corresponding  measurement  equations  are  functionally  independent  of 
each  other.  This  fact  lets  us  reduce  the  dimension  of  the  IS  and  construct  two  4-dim 
IS  Iy(t)  and  Iz{t)  instead  of  one  6-dim  information  set  I(t).  Information  sets  Iy(t) 
and  Iz{t)  are  constructed  separately  for  subspaces  {x,Vr,y,ky}  and  {x,  Vt,  z,  kz}. 

In  order  to  facilitate  the  numerical  operations  with  IS,  it  is  proposed  to  repre¬ 
sent  IS  with  the  help  of  so-called  Grid-Polygon  Approximation  (GPA)  (Fig. 2).  For 
example,  IS  Iy  is  represented  by  a  set  of  its  sections  for  discrete  values  of  x  and  Vt- 

Pij(t)  C  {y,  ky},  i  =  0,...,Nx-l,  j  =  0, . . . ,  Nv  -  1, 

=  4(4  H  {(x,VT,y,ky)  :  x  =  Xi(t),  VT  =  Vj(t)}  , 

where  (xi(t),Vj(t))  is  a  node  of  the  2-dim  grid  on  the  plane  (x,  Vt)-  All  sections 
Pij(t)  are  2-dim  polygons  and  each  polygon  corresponds  to  a  certain  node  of  the 
grid.  The  IS  Iz(t )  is  constructed  in  the  same  way. 

During  the  guidance  the  evolution  of  an  IS  described  by  GPA  consists  of  two 
processes: 

•  motion  of  the  grid’s  nodes; 


•  evolution  of  the  polygonal  sections; 

The  evolution  of  each  polygon  after  getting  a  measurement  at  some  instant  ti  is 
described  by  the  equations: 

where  Hf ,  Hf  are  sets  of  uncertainty  corresponding  to  the  obtained  measurement 
Hi  =  {(y,  ky)  :  —  £  [(u  —  c,  (i;  +  c]}, 

X  i 

H^  =  {(z,  kz)  :  —  G  [(2 1  —  c,  (21  +  c]}. 

xt 

Here  and  later  a  sign  after  time  argument  means  left-side  limit,  for  example 
P-3{ti~)  =  YimP^(t),  ffti. 
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Figure  2:  Grid-Polygon  Approximation  of  an  information  set. 


The  evolution  of  the  information  sets  on  the  interval  [h,  h+i[  is  given  by  the  equations 

Iy{tl+1  )  —  £y  (  ^1/ ),  4(h+ 1  )  —  hb  (  ^2/ )  4  (h  )  , 

where  Cy{-)  and  Cz{-)  -  one-step  extrapolation  operators,  which  can  be  found  from 
(l)-(6)  by  integrating  over  [h,h  +  A]  with  some  controls  u n  and  112V- 


14  = 


—  +  qtA 
Vt 


x'  =  x - —  ln(l  +  o  \i  l  \ / i  /  i A )  -  —  ln(l  +  cyr  14  A),  (9) 

^  Gtrj^ 

IJ  =  y  -  ~  ln(l  +  orVrA)  -  iluA  -  \y{ti)A ,  (10) 

z  =  ~  -  ~ln(l  +  otVtA)  -  u2\A  -  14(h) A,  (11) 

ky  =  ky ,  E  =  k~.  (12) 

Here,  the  sign  ”  1  ”  marks  variables  corresponding  to  the  information  sets  4(h+i  — ) 
and  A(h+i-),  i.e. 

(4,  \'T,  i/,ky)  G  Iy{ti+ 1-)  and  (x,Vr,y,ky)  G  4(h)- 
Analogously, 

(4, 1-4,  4,  E)  G  4(h+i-)  and  (.r,  14,  4)  G  4(h). 
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The  values  Vy(ti)  and  Vz{ti)  are  calculated  by  (5). 

In  this  terms  we  can  define  IS  P  (tn  —  )  extrapolated  to  the  end  of  the  segment 
[t  1 1  t n  [ 

Tly(tn~)  =  £y(uln- 1)  . . .  Cy{un)Iy{ti).  (13) 

The  set  Pz(tn  —  )  is  defined  analogously.  It  is  necessary  to  emphasize  that  extrapo¬ 
lated  IS  P  (tn  —  )  and  Pz(tn  —  )  are  calculated  without  regard  to  future  measurements, 
but  with  regard  to  future  controls  {itu, .  .  .  and  {rt2/,  . .  .  rt2n_i}.  Denote  by 

Pz(tn )  and  Pz(tn)  these  sets  with  regard  to  terminal  controls  U\n  and  u2n- 

4  Auxiliary  impulse  control  problem 

The  aim  of  IS  control  is  to  minimize  expected  terminal  miss  /i  =  (7 Ty,TTz).  This 
terminal  miss  can  be  approximated  by  the  following  equation: 

h2  =  <+^2  =  y2on+^2(n 

where  the  terminal  moment  T  is  defined  by  the  condition 
x(T)  =  0. 

Of  course,  the  accurate  values  of  the  terminal  miss  as  well  as  of  the  moment  T  are 
unknown.  The  idea  is  to  introduce  for  each  point  A  of  the  current  IS  a  miss  function 
tf(A)  -  a  vector- function  with  components  Aly (A)  and  Alz{A).  The  component 
Aly{A)  is  equal  to  the  terminal  miss  value  (by  Y  axis)  obtained  on  controlled  motion 
of  the  system  starting  from  the  point  A  at  instant  t/  and  averaged  over  extrapolated 
IS  Py(T).  The  component  II [(A)  is  defined  analogously.  The  function  rf(A)  is 
taken  as  a  measure  of  the  real  terminal  miss. 

Applying  this  idea  to  the  GPA  of  the  IS  Iy(to)  and  T(t0),  let’s  calculate  first  the 
expected  period  of  the  motion  r°-  for  each  node  of  the  grid  (and  for  each  polygons). 
By  integrating  the  corresponding  dynamics  equations,  we  obtain 

Vr(t)  =  - ,  VM(t)  =  - 

—  +  ar(t  —  to)  w  +  «m(1  —  to) 

Vj  Wo 

To  integrate  equations  (2)  and  (3)  we  calculate  an  auxiliary  integral 
A  1  1  /  \ 

/  -j - ds  —  —  In ( 1  T  olv (t  —  t0) )  (14) 

Jto  -L  |  /  ,  \  OL 

-  +  cds  -  to) 

V 

By  integrating  (2)  from  the  node  (ay,  Vj)  =  (ay(t0),  Vj(t0)),  in  view  of  (14),  we  have 

x(t)  =  xt - —  In [l  +  aMW0(t  -  t0))  -  —  ln(l  +  aTVj (t  -  t0))  •  (15) 

CXm  '  aT  ' 
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(16) 


It  is  clear  that  r°-  is  a  root  of  Eq.  (15)  when  x(t)  =  0,  i.e. 

Xi - —  ln(l  +  (XmWoT?-)  -  —  ln(l  +  aTV3r b)  =  0. 

a*M  v  CKX  v  7 

Assuming  that  the  values  r°-  are  found  by  numerical  solution  of  the  Eq.(16)  for  all 

nodes,  we  can  calculate  the  vector  miss  (TTy}TTz)  for  arbitrary  instant  t  =  ti.  Indeed, 

for  the  instant  ti 

Tij  =$-(!,- to).  (17) 

It  is  assumed  that  if  rj-  <  0  for  some  node  this  node  will  be  excluded  from  the 
further  calculation.  Let  (ay,  V),  y,  ky)  be  a  point  of  the  IS  Iy(ti).  Then  the  miss  nl 
is  calculated  as  follows: 

7T‘y  =  y  -  —  ln(l  +  aTVjT-j)  -  Uik^-J  -  (tk  -  L))  -  Vy{ti)Tll3 , 

°T  tke[tuel] 

where  6l  is  a  lower-bound  estimate  of  the  time  of  control  process  termination  defined 
by  condition  (7).  By  definition, 

el  =  tt  +  A  ^  ,  (is) 

where  [a]  means  an  integer  part  of  a,  A  is  a  time  interval  between  the  observations, 
and  t[  is  a  root  of  the  following  equation 


xLin  -  e - —  ln(l  +  aMVM(ti)4)  -  —  ln(l  +  arV^T^j  =  0, 

(Xm  V  7  CKX  V  7 


where 


xLn  =  V^ax  =  max  Vj(ti). 


Then  the  expression  for  the  averaged  miss  Ilj,  (analogously  n()  takes  the  form 

n(,  =  /  ttL  dx  dVj  dy  dkv,  (19) 

"  QyJlym  V  K  ! 

where  by  a  sign  tilde  we  marked  variables  extrapolated  to  the  moment  6l .  In  par¬ 
ticular  Iy{61)  is  the  IS  Iy(ti)  extrapolated  to  the  moment  6l  and  Qy  is  a  volume  of 
this  extrapolated  IS. 

Now  we  can  formulate  an  auxiliary  control  problem. 

Problem  1.  Let  ti  be  a  current  instant  and  Iy(ti)  be  information  set  after  processing 
measurement  (u.  Let  6l  be  predictable  instant  of  control  process  termination  given 
by  (18).  It  is  required  to  choose  the  impulse  control  u\(tk)  =  Uik,  tk  G  [L,  6l]  in 
order  to  minimize  the  performance  index 

k(“)=  e  <  <2°) 

peEA] 

and  satisfy  the  terminal  condition 

n‘=o.  (2i) 

Problem  2  for  control  process  u2(tk)  =  u2k  is  formulated  analogously. 
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5  Solution  of  auxiliary  control  problem 


It  follows  from  (8)- ( 1 2)  that  the  information  set  extrapolation  (13)  can  be  described 
by  the  following  equations 

1 


Vt  = 


x  = 


1 

—  +  ajT 

Vt 

x - —  ln(l  +  aMVM(ti)T1^  -  —  ln(l  +  aTVTT '), 


y  =  y  -  —  ln(l  +  aTVTT1^  -  uik{rl  -  (tk  -  *;))  -  Vy{ti)r\ 

°T  tke[ti,ol] 

h  —  h 

r^y  —  r^y  •> 


where 


(x,  Vt ,y,ky)  G  Iy(ti ),  {x,  Vt ,y,ky )  G  Iy{6  )• 

Due  to  the  introduced  GPA,  x  =  ay  and  Vt  =  Vj  for  some  node  with  number 
so  t1  equals  to  t1-  from  Eq.  (17).  These  formulas  define  a  change  of  variables  in  the 
integral  (19).  Jacobian  of  this  exchanging  is  equal  to 

Di  =  (l  +  q,tVtt ^ 

In  this  terms  condition  (21)  takes  the  form 

/  [y  “  ln(x  +  aTV3Tlj)  -  Vyitiylj- 

lyVl) 


-  “  (G  -  t))]  (l  +  ocTVjT-^  dx  dVr  dy  dky  =  0.  (22) 

tke[ti,el] 


Here,  the  integral  over  polygons  can  be  calculated  using  Green’s  formula,  the  other 
two  integrations  (over  dx  and  dVr)  can  be  made  using  grid  approximation  and 
Euler’s  formula. 

To  do  this  we  calculate  two  auxiliary  double  integrals  and  J2m  over  two- 
dimensional  polygon  crm,  which  defined  by  a  set  of  its  vertex  {£p,  r]p },  p  =  1, .  .  .  ,  to. 
These  integrals  have  the  forms 

Jm  —  /  J{C-y)dCdih  J2m  =  j  j  C  d^  dy,  C  =  const.  (23) 

&  m  &  m 

We  have  by  Green’s  formula 

Jm  =  J  f(Z-y)d£dy=JZy  (<G  +  dy), 

arn  r  m 
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(24) 


where  Ym  is  a  boundary  of  am.  The  last  integral  is  a  sum  of  m  integrals  J^(p) 

m  m  r 

jl  =  '£Ji(p)  =  !:/{>)(<-({+*(), 

_ _ i  _ _ i  J 


P= i 


7p 


P= 1 

where  each  J^G9)  is  an  integral  over  a  segment  7P  with  ends  {£p,  rjp }  and  {£p+i,  yp+i}, 
p  =  1, .  .  .  ,  to.  It  is  assumed  that  {£m+i,  pm+i}  =  {£i,  71}.  A  segment  7P  admits  a 
parametric  description 

£  =  H^p+i  —  £p)  +  £p5  v  =  Hvp+i  —  Vp)  +  Vp7  A  ^  [o>  i]-  (25) 

Denoting  A^  =  £p+i  —  and  App  =  r/p+1  —  yp,  we  get  with  the  help  of  (25) 

Jm{p)  =  J  £h  (^£  +  dp)  =  (AA 7,  +  £p)(\Anp  +  pp){A^p  +  Aw)  c/A  = 


7p 


>  ,  ^pPp  +  App^p  A^PAPP 

^pVp  +  - o -  +  - o - 


—  +  Ar]P) 


Analogously,  for  the  integral  we  have  by  Green’s  formula 


4  =ffc  4  dr,  =  | 

dm  r 


-vdO  = 

p= i 


where 


(Jr  (J 

•Pnip)  =  77  dr,  —  r,  d£)  =  —  ^PA,JP  —  r,pA,p j. 


(26) 


(27) 


(28) 


7p 


Now  we  calculate  the  integral  (22).  For  this  purpose  we  introduce  the  next 
notations: 


(l  +  aTVjJ 


2  ’ 


=  ^ln(X  +  aT^'ri)’ 


fld=  j  j  {y-kyj 


p?Au) 


(l  +  aTV3TlJ 


dy  dky , 


^■  = 


/  /  K(*0r« 


GAG 


(l  +  aTV3TlJ 


dy  dky , 


where  PT(t;)  is  a  polygon  with  vertices  {yp,  £yp},  p  =  1, .  .  . ,  to,  in  the  plane  (y,  ky 
corresponds  to  the  grid  node  with  number  (i,j)  (Fig. 2). 
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Under  fixed  /,  i,j  we  put 
i  =  V  and  r]  =  kyb\3. 

Then  a  polygon  PF(ti)  transforms  to  a  polygon  am  with  vertices  (^p,gp)  where 
Cp  =  Up  and  ilp  =  kypb *•  •,  p  =  1, .  . .  ,  to,.  In  these  notations 


where  from  (23)  and  it  is  calculated  by  formulas  (24),  (26);  from  (27)  with 
constant  (7  =  1. 

Next,  denote 

d3(k)=  J  J  altJ(k)  dy  dky  =  J2m, 

where  J ^  from  (27)  with  constant  C  =  a\j{k)  and  parameters  £p  =  yp,  r]p  =  kyp. 
The  value  of  this  integral  is  calculated  by  (27)  and  (28). 

As  already  mentioned  above,  we  calculate  the  integral  (22)  by  the  variables  x 
and  Vt  with  the  help  of  grid  approximation  and  Euler’s  formula.  This  gives  in  view 
of  introduced  notations 

Fl-  £  ulkGl(k)  =  0,  (30) 

tke[ti,el] 


where 


F‘  -  E(Ej  +  Fi+ 1  j  +  Fij+1  +  ^i+1  j+l) ’  F‘j  -  fij  - 


hJ 


Gl(k)  —  ^2{g\y{k)  +glt+lj(k)  +  g\j+i(k)  +  g\ 

+1  j+1  (k) 


TO 


Thus,  we  get  a  quadratic  programming  problem  with  the  payoff  function  (20)  and 
linear  condition  (30).  To  solve  this  problem,  let’s  introduce  the  following  Lagrangian 
function 


L(u,  A)  =  J2  <  +  ^\Fl-  £  uikGl(k)  , 
tke[ti,el]  \  tke[ti,el]  J 

where  A  >  0  is  a  Lagrange  multiplier.  It  follows  from  Kuhn- Tucker  theorem  that 
the  optimality  condition  for  the  control  is 


dL(u ,  A) 
duxk 


k  '■  tk  G  [7,  01]. 


(31) 
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From  (31)  and  (30)  we  obtain  the  optimal  control  u*lk  and  parameter  A: 

Gl(k)Fl  x  2  Fl 

uik  Sl  7  A  Sl  ’ 

where 

S‘=  £  [G‘(k)]\ 

peFA] 

The  optimal  control  impulse  at  instant  A  is  given  as 

*  _  G\l)Fl 


6  The  guidance  algorithm  against  a  maneuvering 
target  based  on  the  IS  concept 

So  far  we  considered  only  case  of  target’s  motion  without  any  maneuvers.  In  this 
section,  we  will  consider  a  maneuvering  decelerating  target  lateral  acceleration  which 
can  change  instantly  remaining  within  given  boundaries.  We  will  propose  an  ap¬ 
proach  to  the  treatment  of  possible  target’s  maneuvers  in  the  framework  of  the 
information-set  concept. 


6.1  Simple  one-dimensional  example 

Let’s  consider  one-dimensional  motion  with  model 
x(t)  =  x0  +  a(t), 

where  x0  -  an  initial  position,  a(-)  -  a  term  describing  possible  target’s  maneuver: 
\a(t)\  <  Am,  t  >  t0. 

The  initial  position  x0  is  unknown  but  assumed  to  belong  to  given  interval  I0  which 
defines  an  initial  IS.  The  model  of  measurements  is  very  simple: 

(i  =  x(ti)  +  rji, 

where  T  are  discrete  instants  of  time,  T  =  t0  +  *A,  r/i  is  a  measurement  noise  assumed 
to  be  bound  by  the  following  inequality 

W\  < c- 

After  each  measurement  the  information  set  I~  transforms  as  follows 

u  =  17  n  {* :  ^  g  ic*  -  a  o + c]> 

The  proposed  treatment  of  possible  target’s  maneuver  is  illustrated  by  Fig. 3. 
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Figure  3:  The  behaviour  of  one-dimensional  IS  during  a  target’s  maneuver  and 
proposed  approach  for  its  treatment. 

The  basic  idea  is  to  artificially  increase  size  of  the  IS  during  one-step  exrapolation 
between  two  neighboring  measurements  in  order  to  prevent  possible  IS  degeneration 
due  to  the  maneuver.  In  our  simple  example,  corresponding  equation  of  one-step 
extrapolation  has  a  form 

Ii+ 1  =  [la  —  AmA,  / 1,  +  AmA] , 

where  /&  and  Ia  are  upper  and  lower  ends  of  the  IS  respectively.  Actually,  this 
equation  means  that  IS  is  to  be  inflated  with  maximal  possible  speed  {Am)  on  the 
intervals  between  measurements. 

Let’s  apply  this  approach  to  the  guidance  problems  considered  in  previous  sec¬ 
tions  taking  possible  target’s  maneuvers  into  account. 

6.2  Problem  of  guidance  to  a  maneuvering  target 

The  statement  of  guidance  problem  which  we  will  consider  is  almost  the  same  as  in 
section  2  except  for  modified  equations  of  ’’target-interceptor”  motion. 

In  order  to  describe  target’s  T  maneuvers  it  is  assumed  that  the  target  can 
slightly  change  attack  and  sliding  angles  trying  to  avoid  collision  with  the  inter¬ 
ceptor.  Mathematically  it  may  be  formalized  in  the  next  form.  Suppose  that  in 
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(32) 


equations  (3), (4)  of  ” target-interceptor”  motion 

ky(t )  =  sin^T  /  ay(s)ds , 

Jto 

kz(t )  =  cos  ip  sin  ip  T  /  az(s)ds. 

Jtn 


(33) 


Terms  ay(t),  az(t )  account  for  possible  target’s  maneuvers.  The  only  information  on 
these  functions  is  that  they  are  bounded  by  given  inequalities 


ujV)I  <  -U,  l«V)l  <  4,. 


(34) 


All  the  IS  have  the  same  construction  as  above,  but  their  evolution  on  the  interval 
|_i[  is  given  by  the  equations 


Iy{tl+ 1  )  —  •M.yC'y{ull)Iy{tl)i 

Iz(tl+ 1  )  —  Ai.  z£z  {^21 )  Iz  (tl ) ) 

where  Cy(-)  and  Cy(-)  -  one-step  extrapolation  operators,  which  can  be  found  as 
above  from  (8)- ( 1 2) ,  assuming  that  there  is  no  target’s  maneuver  in  [t/,  t;  +  A[.  Oper¬ 
ators  A4y,  Aiz  perform  linear  transformation  of  the  GPA’s  polygons  in  order  to  take 
target’s  maneuvers  into  account.  As  follows  from  the  model  of  maneuver  (32) — (34) , 
each  point  (y,  ky )  of  a  polygon  is  transformated  into  an  interval  [(y_,  (y+,  k+ )] 

due  to  the  possible  maneuver: 


y  =  y  -  AVL{ vT), 
y+  =  y  +  AyL(Vr ), 


ky  —  ky  Ay  A, 

=ky+  A,,  A, 


(35) 


where  the  function  L(-)  can  be  obtained  by  integrating  equation  (3): 

L(Vt )  = - ttt~  In  (1  +  cdtVtA)  . 

ay  dj  Vt 

Geometrically,  such  transformation  is  nothing  but  stretching  of  the  GPA’s  polygons 
along  a  straight  line  with  slope  L(Vt) / A  as  shown  in  Fig. 4.  The  extrapolated  IS  is 
calculated  in  the  same  way  as  for  non-maneuvering  target,  i.e.  regardless  to  future 
maneuvers  of  the  target  (see  Eq.  (13)).  It  means  that  we  can  use  the  same  statement 
of  the  auxiliary  impulse  control  problem  for  both  maneuvering  and  non-maneuvering 
targets.  Of  course,  the  solution  of  the  problem  will  be  also  the  same. 
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7  Monte  Carlo  test  of  the  algorithm  of  guidance 
to  a  non-maneuvering  target 

In  this  section,  we  present  and  discuss  results  of  Monte  Carlo  tests  which  have  been 
conducted  for  the  algorithm  of  guidance  to  a  non-maneuvering  decelerating  target. 
These  results  regarding  to  accuracy  of  guidance  and  expense  of  controlling  impulses 
will  be  compared  with  the  optimal  solution  of  guidance  problem  with  complete 
information. 

Also  we  will  consider  aspects  of  the  algorithm’s  program  implementation  and 
its  computational  properties,  in  particular,  required  computing  time  and  numerical 
stability  of  the  algorithm. 

7.1  A  priori  data 

The  flow  chart  of  the  guidance  process  simulation  in  the  case  of  constant  lateral 
acceleration  is  shown  in  Fig. 5. 

The  algorithm  was  implemented  in  G  program  module  and  tested  on  PC  Pentium 

III  500  MHz. 

The  algorithm  was  tested  on  a  model  with  parameters  from  Table.  1. 

The  GPA  consisted  of  51x51=2551  nodes: 

Xi  =  •'•()  +  *Ar,  l  /  =  I,,  +  j Ay,  /•  j  0 . 50, 
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Figure  5:  Simulation  of  the  guidance  process  with  constant  lateral  acceleration. 


where 


x0  =  65000  m,  Ax  =  200  m,  Vo  =  4000  m/s,  Ay  =  200  m/s 


The  initial  GPA  polygons  have  been  assumed  to  be  rectangles  defined  by  the  fol¬ 
lowing  inequalities 


Vo  +  tyny  -  Ay 

< 

y 

< 

Vo  +  tiyny  +  Ay, 

~o  +  $zlh  -  A~ 

< 

z 

< 

~o  4“  fi-V-  T  Az, 

sin(<y?  -  Av) 

< 

ky 

< 

sin(<y5  +  Av), 

cos(<^>  +  Av)  sin(<^  -  Av) 

< 

kz 

< 

cos(<y?  —  A^)  sini 

Parameter 

Value 

Parameter 

Value 

Vm 

3000  m/s 

Q'M 

1.5T0-6 

\'t 

4750  m/s 

Q'y 

2.0T0-5 

A 

0.1  s 

V 

-20  grad. 

4' 

20  grad 

Table  1:  Parameters  used  in  the  MG  test. 
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Test  number 

cy,  rad 

cz,  rad 

Number  of  MC  runs 

1 

0.001 

0.001 

1000 

2 

0.002 

0.002 

1000 

3 

0.003 

0.003 

1000 

Table  2:  Values  of  cy,  cz  parameters. 


where  j/q,  Zq  are  initial  positions  which  provide  zero  terminal  miss  on  free  (i.e.  with 
zero  controls)  motion  of  the  ’’target-interceptor”  system,  r]y,  r/z  -  random  values 
uniformly  distributed  over  [—1,  1]  interval.  The  other  parameters  are  the  following 

Sy  =  100  m,  Sz  =  100  m,  Ay  =  3000  m,  A^  =  3000  m, 

Av  =  5.0  grad,  Ay  =  5.0  grad 

The  MC  tests  have  been  conducted  for  three  different  values  of  restrictions  im¬ 
posed  on  measurement  errors  represented  in  Table. 2. 

7.2  Evaluation  of  the  guidance  algorithm  performance 

To  estimate  performance  of  the  guidance  algorithm,  the  following  values  have  been 
calculated  after  each  MC  run: 

1.  Terminal  misses  along  Y  and  Z  axis  My,  Mz. 

2.  Consumption  of  control  impulses,  Ry  and  Rz 

Ry  —  ^  Rz  —  ^ 

f  f 


where  tj  is  the  instant  of  guidance  process  termination. 

The  corresponding  distributions  are  shown  in  Fig. 6  and  Fig. 7.  Filled  histograms 
in  these  figures  show  results  for  the  optimal  solution  of  guidance  problem  with  exact 
information.  Term  ’’exact  information”  means  that  the  initial  parameters  of  the 
target’s  motion  are  taken  from  Monte  Carlo  truth.  After  that  the  deterministic 
problem  of  optimal  control  is  solved.  In  this  problem,  the  performance  index  is  a 
minimization  of  control  expenses  while  terminal  conditions  provide  zero  terminal 
miss. 


7.3  Computational  properties  of  the  algorithm 

In  order  to  understand  the  computational  properties  of  the  proposed  algorithm  we 
have  investigated 
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•  Possibility  of  GPA  degeneration  when  it  prematurely  collapses  and  there  is  no 
non-empty  nodes  at  the  end  of  the  guidance. 

•  Convergence  of  GPA  to  the  true  values  of  x0  and  Vt  parameters.  In  case 
of  convergence  all  non-empty  nodes  of  GPA  should  be  located  around  point 
corresponding  to  the  true  values  of  x0  and  Vt- 

•  Computational  speed  of  the  algorithm  and  timing  of  different  parts  of  the 
algorithmic  procedure. 

The  histograms  for  the  mean  computing  time  and  number  of  non-empty  GPA 
nodes  at  the  end  of  the  guidance  are  presented  in  Fig. 8  -  10. 

One  can  see  that  zero  bin  of  the  histograms  for  the  final  number  of  GPA  nodes 
is  empty.  It  means  that  during  all  3000  (3x1000)  MC  runs  there  were  no  cases  of 
premature  collapse  of  GPA  regardless  to  the  measurement  accuracy. 

The  mean  computing  time  and  partial  timing  are  summarized  in  Table. 3. 


Values  of  cy,  cz,  rad 

0.001 

0.002 

0.003 

Total  mean  time,  s 

4.06 

4.76 

5.26 

Information  set  filter,  s 

1.19 

1.33 

1.43 

Calculation  of  control,  s 

2.09 

2.46 

2.76 

One-step  GPA  extrapolation,  s 

0.75 

0.94 

1.04 

Table  3:  MC  test  results  on  computational  speed  of  the  algorithm 

The  correct  convergence  of  the  GPA  is  illustrated  by  Fig.  11  -  13.  On  these  plots, 
the  average  GPA  profiles  are  shown  for  different  instants  of  time.  By  definition,  the 
GPA  profile  is  a  two-dimensional  distribution  calculated  as  follows 

Cruris 

E  Ek(i,j,t) 

k= 0 


where  Nruns  -  the  number  of  MC  runs  taken  for  averaging,  Ek(i,j,t )  equals  to  1 
if  GPA  node  (i,j)  at  instant  t  is  non-empty  during  MC  run  number  k,  and  equals 
to  0  if  this  node  is  empty.  The  true  values  of  x0}  Vt  parameters  are  marked  by 
haircrosses. 


Vi-jV) 
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Terminal  miss  My,  m 


Control  expense  Ry,  m/s 


Terminal  miss  My,  m 


Control  expense  Ry,  m/s 

ure  6:  MG  test  results  for  cy 
and  cy  =  c~  =  0.002  rad. 


Terminal  miss  Mz,  m 


Control  expense  Rz,  m/s 


Terminal  miss  Mz,  m 


Control  expense  Rz,  m/s 

c~  =  0.001  rad  (upper  two  plots) 
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Figure  7:  MC  test  results  for  cy  =  c~  =  0.003  rad. 


Computing  time,  s 


Final  number  of  non-empty  GPA  nodes 

Figure  8:  MC  test  results  for  cy  =  c~  =  0.001  rad. 
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Average  G P A  profile,  t  =  7  s 
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Average  GPA  profile,  t=9  s 


Average  GPA  profile,  t  =  1 1  s 


Final  rjyerage  GFA  profile 


Figure  11:  Average  GPA  profiles  for  cy  =  c~  =  0.001  rad. 
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Average  GPA  profile,  t  =  1  1  s 


Final  average  GPA  profile 
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Figure  12:  Average  GPA  profiles  for  cy  =  c~ 


0.002  rad. 
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Figure  13:  Average  GPA  profiles  for  cy  =  c~  =  0.003  rad. 
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7.4  Resume 

The  results  obtained  during  MC  test  of  the  proposed  guidance  algorithm  show  that 


•  The  algorithm  provides  precise  and  unbiased  guidance,  the  distributions  of 
terminal  misses  are  symmetrical  without  pronounced  tails  for  both,  OY  and 
OZ,  control  directions. 

•  The  consumption  of  control  impulses  during  guidance  depends  on  given  mea¬ 
surement  errors  and  tends  to  the  optimal  solution  with  complete  information 
when  the  measurement  precision  is  increased. 

•  The  proposed  algorithm  is  numerically  stable,  no  cases  of  GPA  degeneration 
have  been  observed  during  3000  MC  runs  with  different  values  of  measurement 
errors. 

•  The  final  average  distribution  of  non-empty  GPA  nodes  has  pronounced  single 
peak  corresponding  to  the  true  values  of  x0  Vt  target  motion  parameters.  It 
means  that  the  algorithm  can  identify  these  unknown  parameters  correctly 
and  GPA  finally  converges  to  them.  The  speed  of  such  convergence  is  higher 
for  more  accurate  measurements. 

•  The  mean  computing  time  of  the  algorithm  highly  depends  on  the  measure¬ 
ment  accuracy  and  decreases  when  the  accuracy  is  increased.  It  can  be  ex¬ 
plained  by  fast  convergence  of  GPA  in  case  of  more  accurate  measurements. 
In  turn,  low  number  of  remaining  GPA  nodes  leads  to  the  faster  work  of  the 
algorithm. 
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8  Monte  Carlo  test  of  the  algorithm  of  guidance 
to  a  maneuvering  target 

In  this  section  we  present  and  discuss  results  of  Monte  Carlo  tests  which  have  been 
conducted  for  the  algorithm  of  guidance  to  a  maneuvering  decelerating  target. 

As  we  mentioned  above,  it  is  assumed  that  the  target  can  perform  evasive  ma¬ 
neuvers  by  changing  its  sliding  and  attack  angles.  In  turn,  it  leads  to  additional 
non-zero  lateral  acceleration  which  change  the  direction  of  the  target’s  velocity  vec¬ 
tor  in  chosen  coordinate  system.  In  order  to  describe  this  mechanism  of  the  target’s 
maneuver  in  the  framework  of  IS-based  approach  we  assume  that  lateral  acceleration 
of  the  target  can  change  instantly  remaining  within  given  boundaries. 

At  the  beginning,  we  confirm  by  Monte  Carlo  simulation  that  the  algorithm  of 
guidance  to  a  non-maneuvering  target  considered  above  fails  in  the  case  when  the 
targets  makes  an  evasive  maneuver.  As  was  shown,  not  taking  a  maneuver  into 
account  causes  fast  degeneration  of  CPA’s  nodes  and  large  terminal  misses  which 
allows  the  target  to  avoid  an  intercept. 

Then  we  show  that  the  algorithm  proposed  in  section  6  provides  reliable  guidance 
to  a  maneuvering  target.  The  distributions  of  terminal  misses  for  this  algorithm  are 
significantly  better  than  that  for  the  algorithm  of  guidance  to  a  non-maneuvering 
target.  Results  regarding  to  expense  of  controlling  impulses  are  compared  with  the 
optimal  solution  of  guidance  problem  with  complete  information.  Also  we  consider 
aspects  of  the  algorithm’s  computational  properties,  in  particular,  required  comput¬ 
ing  time  and  numerical  stability  of  the  algorithm. 

8.1  A  priori  data 

The  flow  chart  of  the  guidance  process  simulation  in  the  case  of  maneuvering  target 
is  shown  in  Fig.  14.  In  comparison  with  the  chart  in  Fig. 5,  this  one  has  additional 
blocks  ’’Simulation  of  target  maneuver”  and  ’’Linear  transformation  of  CPA  nodes”. 

This  algorithm  was  implemented  in  C  program  module  and  tested  on  PC  Pentium 
III  500  MHz  too. 

All  the  a  priori  data  were  the  same  as  in  the  case  of  non-maneuvering  target. 
Namely,  the  CPA  consisted  of  51x51=2551  nodes: 

Xi  =  x0  +  iAx,  VTj  =  Co  +  j Ay,  i,j  =  0, . . . ,  50, 

where 

x0  =  65000  m,  A^  =  200  m,  Vo  =  4000  m/s,  Ay  =  200  m/s 
The  initial  CPA  polygons  have  been  assumed  to  be  rectangles  defined  by  the  fol- 
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lowing  inequalities 

yo  +  SyVy-Ay  <  y  <  yo  +  8yrjy  +  Ay, 

~o  +  S~i]z~A~  <  z  <  z*  +  5zi]~  +  A~, 

sin(<y5  —  A^)  <  ky  <  sin(<y5  +  A^), 

cos(<^>  +  Av)  sin(<^>  -  Av)  <  k~  <  cos(<^>  -  Av)  sin(<^>  +  Av), 

where  y£,  ~q  are  initial  positions  which  provide  zero  terminal  miss  on  free  (i.e.  with 
zero  controls)  motion  of  the  ’’target-interceptor”  system,  i]y,  y~  -  random  values 
uniformly  distributed  over  [—1,  1]  interval.  The  other  parameters  are  the  following 

8y  =  100  m,  8z  =  100  m,  Ay  =  3000  m,  A~  =  3000  m, 

Av  =  5.0  grad,  Ay  =  5.0  grad 

The  algorithm  was  tested  on  a  model  with  parameters  from  Table.  1  under  three 
different  values  of  restrictions  imposed  on  measurement  errors  represented  in  Table.  2. 

As  was  mentioned  above,  we  consider  the  following  mechanism  of  the  target’s 
maneuver.  At  som#  instants  iy,  the  angles  0,  Lp  which  del  ermine  a  position  of 
the  target’s  velocity  vector  start  changing  (Fig. 15).  The  rates  0,  Lp  are  assumed  to 
be  constant.  The  maneuver  durations  are  Ty,  70  for  angles  0,  Lp  correspondingly. 

For  the  Monte  Carlo  simulation  described  here,  the  following  values  have  been 
taken  (Table  4). 
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Figure  15:  Evolution  of  0  and  p>  during  target’s  maneuver. 


Parameter 

Value 

Parameter 

Value 

t'lp 

5.0  s 

t(p 

5.0  s 

t4, 

0.2  s 

0.2  s 

4' 

0.2  grad/s 

V 

0.2  grad/s 

Table  4:  Maneuver  parameters  taken  for  the  Monte  Carlo  simulation. 


In  order  to  take  possible  target’s  maneuvers  into  account,  the  GPA  polygons  have 
been  subjected  to  a  special  linear  transformation  in  addition  to  the  one-step  extrap¬ 
olation  of  their  vertices  in  accordance  with  (35).  The  transformation  parameters 
A~  have  been  takes  as  (Table  5) 


Parameter 

Value 

Parameter 

Value 

Ay 

3.5T0-4  s"1 

Az 

3.5T0-4  s"1 

Table  5:  Parameters  used  for  transformation  of  GPA’s  polygons. 


8.2  Evaluation  of  the  guidance  algorithm  performance 

To  estimate  performance  of  the  guidance  algorithm,  we  have  used  the  same  method 
as  in  section  7.2,  namely,  the  following  values  have  been  calculated  after  each  MG 
run: 

1.  Terminal  misses  along  Y  and  Z  axis  My,  M~. 

2.  Consumption  of  control  impulses,  Ry  and  R~ 

Ry  —  0  '  Rz  —  0  ' 


where  tj  is  the  instant  of  guidance  process  termination. 
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8.3  Computational  properties  of  the  algorithm 

As  in  the  case  of  non-maneuvering  target,  in  order  to  evaluate  the  computational 
properties  of  the  proposed  algorithm  we  have  investigated 


•  Possibility  of  GPA  degeneration  when  it  prematurely  collapses  and  there  is  no 
non-empty  nodes  at  the  end  of  the  guidance. 

•  Convergence  of  GPA  to  the  true  values  of  x0  and  Vt  parameters.  In  case 
of  convergence  all  non-empty  nodes  of  GPA  should  be  located  around  point 
corresponding  to  the  true  values  of  x0  and  Vt- 

•  Computational  speed  of  the  algorithm  and  timing  of  different  parts  of  the 
algorithmic  procedure. 


The  correct  convergence  of  the  GPA  have  been  checked  by  investigating  so- 
called  average  GPA  profiles  for  different  instants  of  time.  These  two-dimensional 
distributions  are  calculated  as  follows 


'Pp-(t) 


E 


k= 0 


Ekij-i  jV) 


Nr. 


•) 


where  Nruns  -  number  of  MC  runs  taken  for  averaging,  Ek(i,j,t )  equals  to  1  if  GPA 
node  (i,j)  at  the  moment  t  is  non-empty  during  MC  run  number  k,  and  equals  to 
0  if  this  node  is  empty. 


8.4  Results  of  the  Monte  Carlo  simulation 

As  was  mentioned  above,  we  have  investigated  first  the  behaviour  of  the  guidance  al¬ 
gorithm  for  non-maneuvering  target  in  the  case  when  a  target  performs  a  maneuver. 
Obtained  distributions  of  terminal  misses  and  average  GPA  profiles  are  presented  in 
Fig.  16  -  18  for  three  different  values  of  restrictions  imposed  on  measurement  errors 
(Table  2).  The  true  values  of  x0}  Vt  parameters  are  marked  by  haircrosses. 

Histograms  for  the  final  number  of  GPA  nodes  show  that  premature  collapse 
of  GPA  have  happened  very  often  for  the  guidance  algorithm  under  consideration. 
Observed  large  terminal  misses  are  explained  by  complete  degeneration  of  the  CPA’s 
nodes  due  to  the  target’s  maneuver  which  is  confirmed  by  corresponding  average 
GPA  profiles. 

On  the  contrary,  the  algorithm  proposed  for  maneuvering  target  provides  very 
reliable  guidance  against  a  maneuvering  target.  The  distributions  of  terminal  misses 
and  consumption  of  control  impulses  are  presented  in  Fig.  19  -  20.  Filled  histograms 
in  these  figures  show  results  for  the  optimal  solution  of  guidance  problem  with  exact 
information.  As  was  mentioned  above,  term  ’’exact  information”  means  that  the 
initial  parameters  of  the  target’s  motion  are  taken  from  Monte  Carlo  truth.  After 
that  the  deterministic  problem  of  optimal  control  is  solved.  In  this  problem,  the 
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performance  index  is  a  minimization  of  control  expenses  while  terminal  conditions 
provide  zero  terminal  miss. 

The  histograms  for  the  mean  computing  time  and  number  of  non-empty  GPA 
nodes  at  the  end  of  the  guidance  are  presented  in  Fig. 21  -  23. 

One  can  see  that  zero  bin  of  the  histograms  for  the  final  number  of  GPA  nodes 
is  empty.  It  means  that  during  all  3000  (3x1000)  MC  runs  there  were  no  cases  of 
premature  collapse  of  GPA  regardless  to  the  measurement  accuracy. 

The  mean  computing  time  and  partial  timing  are  summarized  in  Table  6. 


Values  of  cy,  cz,  rad 

0.001 

0.002 

0.003 

Total  mean  time,  s 

10.81 

13.35 

16.27 

Polygons  intersections,  s 

2.88 

3.54 

4.31 

Calculation  of  control,  s 

4.42 

5.46 

6.65 

One-step  GPA  extrapolation,  s 

3.48 

4.32 

5.28 

Table  6:  MC  test  results  on  computational  speed  of  the  algorithm 

The  correct  convergence  of  the  GPA  is  illustrated  by  the  average  GPA  profiles 
shown  in  Fig. 24  -  26.  The  true  values  of  x0}  Vt  parameters  are  marked  by  haircrosses. 

8.5  Resume 

The  results  obtained  during  MC  test  of  the  proposed  guidance  algorithm  show  that 

•  The  algorithm  provides  precise  and  unbiased  guidance,  the  distributions  of 
terminal  misses  are  symmetrical  without  pronounced  tails  for  both,  OY  and 
OZ,  control  directions. 

•  The  consumption  of  control  impulses  during  guidance  depends  on  given  mea¬ 
surement  errors  and  tends  to  the  optimal  solution  with  complete  information 
when  the  measurement  precision  is  increased. 

•  The  proposed  algorithm  is  numerically  stable,  no  cases  of  GPA  degeneration 
due  to  the  target’s  maneuver  have  been  observed  during  3000  MC  runs  with 
different  values  of  measurement  errors. 

•  The  final  average  distribution  of  non-empty  GPA  nodes  has  pronounced  single 
peak  corresponding  to  the  true  values  of  x0  Vt  target  motion  parameters.  It 
means  that  the  algorithm  can  indentify  these  unknown  parameters  correctly 
and  GPA  finally  converges  to  them.  The  speed  of  such  convergence  is  higher 
for  more  accurate  measurements. 
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•  The  mean  computing  time  of  the  algorithm  highly  depends  on  the  measure¬ 
ment  accuracy  and  decreases  when  the  accuracy  is  increased.  It  can  be  ex¬ 
plained  by  fast  convergence  of  GPA  in  case  of  more  accurate  measurements. 
In  turn,  low  number  of  remaining  GPA  nodes  leads  to  the  faster  work  of  the 
algorithm. 

9  Conclusion 

The  three-dimensional  problem  of  guidance  with  incomplete  information  has  been 
considered.  To  solve  the  problem  an  approach  concerned  with  a  construction  of 
information  sets  (IS)  has  been  used.  In  terms  of  IS  the  original  problem  has  been 
reduced  to  an  auxiliary  impulse  control  problem  and  the  guidance  algorithms  against 
high-speed  non-maneuvering  and  maneuvering  targets  has  been  developed.  To  fa¬ 
cilitate  numerical  operations  with  the  introduced  IS,  Grid-Polygon  Approximation 
of  the  IS  has  been  proposed.  To  estimate  performance  of  the  guidance  algorithms, 
Monte  Carlo  (MC)  simulation  has  been  accomplished.  The  MC  simulation  has  been 
demonstrated  numerical  stability  of  the  algorithms  and  their  potential  for  real-time 
data  processing. 


10  Outlook 

In  our  opinion,  the  future  improvements  on  the  proposed  algorithms  will  be  mainly 
concerned  with  increasing  their  computational  speed.  Such  an  improvement  can  be 
done  in  the  different  ways. 

First,  we  can  use  more  efficient  numerical  approximation  of  the  introduced  IS. 
For  example,  new  nodes  can  be  added  to  the  GPA  dynamically  in  order  to  keep  their 
total  number  more  or  less  constant  in  the  course  of  the  guidance.  Also  it  might  be 
better  to  describe  the  sections  of  the  IS  by  not  simply  polygons  but  by  combinations 
of  straight-line  segments  and  segments  of  some  curves  of  higher  order. 

Second,  as  we  mentioned  above,  the  MC  tests  show  that  the  number  of  GPA 
nodes  remains  quite  big  during  all  guidance,  slowing  the  algorithms  significantly 
(one  run  takes  about  7-10  sec).  In  turn,  this  number  highly  depends  on  mea¬ 
surement  accuracy  and  decreases  when  the  accuracy  is  increased.  It  is  natural, 
because  the  higher  measurement  accuracy  is,  the  better  observability  the  system 
’’target-interceptor”  has.  Measurement  accuracy  depends  on  inner  parameters  of  IR 
array  seeker  only  and  cannot  be  improved  by  interceptor’s  impulse  control  law.  In 
this  connection  there  is  the  main  problem:  to  reduce  the  number  of  remained  GPA 
nodes  under  given  measurement  accuracy.  To  solve  this  problem  we  need  to  improve 
observability  of  the  system  ’’target-interceptor”. 

It  is  well  known  that  in  problems  of  guidance  with  nonlinear  (bearings/range, 
bearings-only)  measurements,  the  observability  explicitly  depends  on  guidance  tra¬ 
jectory.  So,  there  is  an  idea  to  try  so-called  trajectory  control  over  observations 


33 


[7],  [8]  to  construct  a  guidance  trajectory  of  special  shape  that  makes  ’’target- 
interceptor”  system  more  observable  and  hence  improves  quality  of  TMA  solution 
that  is  of  vital  importance  in  the  case  of  engagement  with  a  maneuvering  target. 
The  key  element  of  this  idea  is  to  add  into  the  performance  index  of  the  optimal 
impulse  control  problem  the  special  term  that  has  the  effect  of  making  the  target’s 
motion  more  observable.  The  performance  index  is  taken  to  be  the  weighted  sum 
of  an  integral  penalty  on  the  magnitude  of  the  interceptor  controls  and  the  trace 
of  the  Fisher  information  matrix  (observability  Gramian).  This  performance  index 
should  be  minimized  under  terminal  conditions  requiring  that  the  average  terminal 
miss  be  equal  to  zero.  The  averaged  miss  is  defined  as  a  ratio  of  the  integral  of  the 
terminal  miss  function  over  the  IS  and  the  IS  phase  volume.  In  our  opinion,  this  is 
an  interesting  problem  for  future  investigations. 
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Figure  16:  MG  test  results  for  cy  =  c~  =  0.001  rad. 
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Figure  17:  MG  test  results  for  cy  =  c~  =  0.002  rad. 
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Figure  18:  MG  test  results  for  cy  =  c~  =  0.003  rad. 
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Figure  19:  MG  test  results  for  cy  =  c3 
0.002  rad. 
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Figure  20:  MC  test  results  for  cy  =  c~  =  0.003  rad. 
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Figure  21:  MC  test  results  for  cy  =  c~  =  0.001  rad. 
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Figure  24:  Average  GPA  profiles  for  cy  =  c~  =  0.001  rad. 
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Figure  25:  Average  GPA  profiles  for  cy  =  c~  =  0.002  rad. 
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Figure  26:  Average  GPA  profiles  for  cy  =  c~  =  0.003  rad. 
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